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ABSTRACT 

T Tauri stars are young, low mass, pre-main sequence stars surrounded by an accretion disk. These objects 
present strong magnetic activity and powerful magnetic reconnection events. Strong shocks are likely associ- 
ated with fast reconnection in the stellar magnetosphere. Such shocks can accelerate particles up to relativistic 
energies. We aim at developing a simple model to calculate the radiation produced by non-thermal relativistic 
particles in the environment of T Tauri stars. We want to establish whether this emission is detectable at high 
energies with the available or forthcoming 7-ray telescopes. We assume that particles (protons and electrons) 
pre-accelerated in reconnection events are accelerated at shocks through Fermi mechanism and we study the 
high-energy emission produced by the dominant radiative processes. We calculate the spectral energy dis- 
tribution of T Tauri stars up to high-energies and we compare the integrated flux obtained with that from a 
specific Fermi source, IFGL J1625.8-2429c, that we tentatively associate with this kind of young stellar ob- 
jects (YSOs). We suggest that under reasonable general conditions nearby T Tauri stars might be detected at 
high energies and be responsible for some unidentified Fermi sources on the Galactic plane. 
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stars: winds, outflows 
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1. INTRODUCTION 

The Fermi Gamma-Ray Space Telescope was launched on 
2008 June 1 1 and is the successor of the Compton Gamma 
Ray Observatory. It has two 7-ray instruments on board, the 
Large Area Telescope (LAT) and the Gamma-ray Burst Mon- 
itor (GBM). The first one, a pair-production telescope, is the 
successor of the Energetic Gamma Ray Experiment Telescope 
(EGRET). LAT has a sensitivity ^ lO"** ph s~' cm"^ in the en- 
ergy range from 20 MeV to over 300 GeV. 

Based on the first eleven months of survey, the Fermi First 
Year Catalog has been published (Abdo et al. 2010). The cat- 
alog consists of 1451 sources with statistical significances of 
4cr or higher This catalog provides potential identifications, 
when available. Nearly 630 Fermi sources remained uniden- 
tified (Abdo et al. 2010). A fraction of these unidentified 
sources lay on the Galactic plane (latitudes < 5°). They may 
be pulsars, supernova remnants (SNR), the effect of strong 
winds from massive stars, high-mass or low-mass X-ray bi- 
naries (HMXBs or LMXBs), or stellar clusters, among other 
possibilities (e.g. Romero, Benaglia, & Torres 1999). 

Munar-Adrover, Paredes, & Romero (2011) have crossed 
the Fermi First Year Catalog with catalogs of known YSOs, 
in order to identify those protostars that might emit 7-rays. 
They conclude that 72% of the candidates obtained by spatial 
correlation should be 7-ray sources with a confidence above 
5(7. Massive YSOs have already been claimed to be 7-ray 
sources (e.g. Araudo et al. 2007, Bosch-Ramon et al 2010). 
However, such a claim has not been done yet for low-mass 
protostars. 
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In this paper we shall argue that these stars, in particular 
T Tauri stars, can be faint but sometimes detectable 7-ray 
sources. We will focus our attention on the physical processes 
that can generate 7-ray emission in T Tauri magnetospheres. 
Specifically, we shall discuss whether these protostars can be 
responsible for sources like IFGL J1625.8-2429c. In the next 
section we present a brief introduction to T Tauri stars and we 
describe the general physical scenario. In Section 3 we pro- 
vide the details of our model. In Section 4 we present our 
calculations applied to a specific group of T Tauri stars and 
the main results. In Section 5, we discuss these results and, 
finally, we close with some brief conclusions. 

2. SCENARIO 

T Tauri stars are low-mass stars (M < 3 Mq) in their early 
stages of evolution. Classical T Tauri stars have typical K-M 
specti-al types with Teff -- 3000-5000 K (Montmerle & Andre 
1989). Values of the stellar radius are ^ 2-3 Rq. As pro- 
genitors of solar-like stars, they are subject of intense study. 
They are usually optically visible and the youngest objects 
drive bipolar outflows. The protostars are born in the collapse 
of molecular clouds. The cloud material has high angular mo- 
mentum and therefore a circumstellar accretion disk is formed 
(e.g. Bekwith et al. 1990). These accretion disks are trun- 
cated in the vicinity of a co-rotation radius where the magnetic 
pressure equals the gas pressure in the disk (see Fig. [U. In- 
frared observations allow to infer inner radii of 0.07-0.54 AU 
(Muzerolle et al. 2003), consistent with the co-rotation radius. 
UV emission has been used to estimate the mass accretion 
rate, yielding typical values of ^ 10"*^ Mq year"^ (Gomez de 
Castro & Lamzin 1999; Johns-Ki-ull et al. 2000). The infall 
velocity can almost reach the free-fall speed. The existence of 
a large amount of dust in a flattened disk is confirmed by in- 
tense emission from the micrometer to the millimeter bands. 
The disks can now be directly imaged (e.g. Dutrey, Guil- 
loteau, & Simon 1994). In addition, these objects drive strong 
winds with mass loss rate M ^ 10~^Mq year"' and velocity 
Vw -- 200 km s"' (Feigelson & Montmerle 1999). 
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X-ray studies indicate particle number densities of the ac- 
creting plasma of about 10'^ cm"-' (Gunther et al. 2007). 

Variable thermal X-ray emission is also detected from T 
Tauri stars in the keV band. Luminosities are found to be in 
the range ^ lO-'' - lO-'-' erg s"'. This emission comes from 
a high density plasma at a typical temperature of ^ 10^ K; 
flares with temperatures ^ 10 K have been detected (Tsuboi 
et al. 1998). These flares have durations of ~ 10^ - lO'* s. 
Such events are expected to occur in magnetic flux tubes with 
spatial extent of ^ 10'" - 10" cm (e.g Hayashi et al. 1996). 
Models for the X-ray activity based on the interaction of the 
stellar object and the circumstellar disk have been proposed 
by several authors (e.g. Hayashi et al. 1996, Birk et al. 2000). 

These X-rays flares are widely considered as upscaled ver- 
sions of solar flares. The rapid heating and cooling of plasma 
and acceleration of particles must arise from efficient MHD 
processes, such as solar-type magnetic reconnection events in 
twisted flux tubes that connect the central object and the cir- 
cumstellar disk (Birk et al. 2000). Magnetic reconnection 
is one of the fundamental processes in astrophysical plasmas 
because it explains large-scale, dynamic releases of magnetic 
energy. It is essentially a topological reconfiguration of the 
magnetic field caused by a change in the connectivity of the 
field lines. It is this change which allows the release of stored 
magnetic energy, which in many situations is the dominant 
source of free energy in a plasma. Several works have been 
done on particle acceleration through magnetic reconnection 
( e.g. Schopper, Lesch, & Birk 1998, Zenitani & Hoshino 
2001, de Gouveia Dal Pino et al. 2010). Strong shocks result- 
ing from supersonic plasma ejections are the likely outcome 
of massive reconnection in T Tauri magnetospheres. Such 
shocks can in principle accelerate particles up to relativistic 
energies through Fermi mechanism. 

The expected values of magnetic field in T Tauri stars are 
^ 1 kG (e.g. Johns-Krull 2007) and the field structure is com- 
plex and multipolar, as in the Sun. For simple flare models, 
quantitative properties of large-scale magnetic field structures 
can be inferred (e.g. Montmerle et al. 1983, Walter & Kuhi 
1984). Assuming equipartition conditions, B^ /^n = ItigkT , 
the magnetic field strength in the magnetosphere is approxi- 
mately 10^ G (Feigelson & Montmerle 1999). 

3. MODEL 

3.1. Particle acceleration 

Diffusive shock acceleration does not work for slow-mode 
shocks, consequently is thought that diffusive shock acceler- 
ation is not important for magnetic reconnection that involves 
slow-mode shocks (e.g. Priest & Forbes 2000). Under some 
conditions, however, any obstacle which obstructs the outflow 
will create a fast-mode shock (Tsuneta & Naito 1998). In a T 
Tauri magnetosphere the obstacles might be clumps from the 
strong stellar wind; considerable observational evidence sup- 
ports the idea that the wind structure is clumpy (e.g. Owocki 
& Cohen 2006). Moreover, shock acceleration is now applied 
to the outflow regions of coronal magnetic reconnection sites, 
where first-order Fermi mechanism at the standing fast shock 
is a leading candidate (Aschwanden 2008). The shocks are 
expected to accelerate charged particles up to high energies 
by a Fermi-like diffusive process (e.g. Drury 1983). 

Additionally, recent extensive 3D numerical simulations 
performed by Kowal et al. (2011) show that within con- 
tracting magnetic islands or current sheets, charged particles 
are accelerated by a first-order Fermi process while outside 
the current sheets and islands the particles experience mostly 



drift acceleration due to magnetic field gradients. These re- 
sults are supported by observations of solar flares that suggest 
that magnetic reconnection should be first slow in order to 
ensure the accumulation of magnetic flux and then suddenly 
become fast to allow a rapid energy release (e.g. Lazarian 
& Vishniac 1999). Particles scattered by turbulence between 
converging magnetic mirrors formed by oppositely-directed 
magnetic fluxes moving towards each other at the fast recon- 
nection speed can undergo diffusive acceleration without the 
requirement of strong shock formation (Kowal et al. 201 1). 

Independently of the details of the acceleration mechanism, 
we assume that a population of non-thermal relativistic parti- 
cles, electrons and protons, is injected into the magnetosphere 
of the T Tauri star. These particles will interact with the large 
scale magnetic fields, with the existing radiation fields, and 
with the magnetospheric plasma, producing non-thermal elec- 
tromagnetic radiation. The size of the acceleration region is 
the spatial scale where reconnection takes place, i.e the flux 
tube length. 

The power available in the magnetized system is 
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where A is the area of the flux tube, of length / = 10" cm and 
aspect ratio '-^ 0.1 Z (Feigelson & Montmerle 1999), and va 
is the Alfven speed, va = ^JB^ /{A-KnipVi), with n the particle 
density and nip the proton mass. For B = 2x 10^ G, we get L 
~ 10^^ erg s~'. We assume that 10% of this power is released 
in the reconnection processes. In turn, a fraction ^rei ^ 10% 
of this power goes to relativistic particles. These values are 
in accordance to inferred values for the Sun, and can be con- 
sidered even conservative. For instance, in large solar flares 
the accelerated particles contain up to 10-50% or more of the 
total energy released, whereas in gradual events ^ 10% of the 
total power goes to the accelerated particles (see Lin 2008 and 
references therein). 

The efficiency of non-thermal acceleration in the magne- 
tized plasma is 
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mimicking the efficiency for standard first order Fermi accel- 
eration theory behind shocks^ (Drury 1983, Vila & Aharonian 
2009); here D is the particle difussion coefficient, r^ is the 
particle gyroradius and Vrec is the reconnection speed. If D 
is in the Bohm limit DBohm = ''gc/3. The reconnection speed 
in violent reconnection events satisfy Vrec ^ va (Lazarian & 
Vishniac 1999, Kowal et al. 2009), so we assume Vrec = 0.6va, 
that gives an efficiency ry ^ lO"''. This efficiency, although not 
very high, will allow maximum energies well into the 7-ray 
domain, as we will see. 

3.2. Energy losses 

The relativistic electrons and protons lose energy through 
different interactions. The maximum energy that they can 
reach depends on the processes of energy loss, the available 

' This assumption liolds only if there is a first-order Fermi mechanism 
occurring within the reconnection zone caused by two converging magnetic 
fluxes of opposite polarity which move to each other with v^^c ■ This mecha- 
nism of first-order Fermi acceleration occurring within a reconnection site 
was first proposed and described analyticafly in de Gouveia Dal Pino & 
Lazarian (2005) and successfully tested through 3D MHD simulations by 
Kowal, de Gouveia Dal Pino, & Lazarian (201 1). 




Fig. 1. — Sketch of a T Tauri star adapted from Feigelson & Montmerle 
(1999). 

space for effective acceleration and on the acceleration rate. 
This rate is t~^^ = -qecB/E for first order Fermi mechanism. 

The main loss mechanism for electrons in the magneto- 
sphere of T Tauri stars is synchrotron radiation. Other rele- 
vant loss processes are relativistic Bremsstrahlung — the lep- 
tons interact with ions of material in the magnetosphere — , 
and inverse Compton scattering (IC) — electrons scatter up 
the ambient photons — . The main ambient photon field is the 
thermal X-ray emission, which corresponds to a temperature 
of ^ 10^ K, and has a typical luminosity (in quiescence) of 
~ 6 X 10^"^ erg s"' according to the observed values (see Table 

El). 

Figure |2] shows the cooling rates and the acceleration rate 
for electrons in the case of the parameters shown in Table 1 . 
These parameters are fixed in accordance to the typical values 
discussed above for T Tauri stars (see e.g. Feigelson & Mont- 
merle 1999). Particles can also escape by convection from the 
acceleration region with a rate t~^^^ = v^/l, where Vw is the 
wind velocity. The maximum energy for electrons is obtained 
equating the cooling rates with the acceleration rate. The re- 
sult is ^ 2 GeV. It can be seen that electrons rapidly cool by 
synchrotron mechanism with timescales of less than 1 s. 

The radiative losses for protons are p^ interactions with 
the X-ray radiation field, synchrotron and pp inelastic col- 
lisions. The latter losses are the most important for these par- 
ticles. Figure [3] shows the cooling rates and the acceleration 
rate for protons. The maximum energy of these particles is 
much greater than the maximum energy of electrons, ^10 
TeV. Protons of such an energy satisfy the Hillas criterion for 
a field of 2 x 10^ G. Acceleration timescales at the higher en- 
ergies are ^ 10^ s, larger than the typical durations of flaring 
X-ray episodes (^ 10"* s). Hence, maximum energies of ^ 10 
TeV are more realistic in a steady state approximation to the 
7-ray active period. 



Fig. 2. — Acceleration and cooling rates for electrons. The identification 
of the different curves is given in the box in the upper right comer. 
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Fig. 3. — Acceleration and cooling rates for protons. 



TABLE 1 
Parameters 



Parameter [units] 




values 


/ Flux tube length [cm] 




10" 


v'A Alfven speed [cm s"' ] 




6x10' 


'q Acceleration efficiency 




io-« 


a Hadron-to-lepton energy 


ratio 


100 


^lei Fractional content of relativistic particles 


10-1 


a Particle injection index 




2.2 


I'w Wind velocity [cm s"'] 




2x10' 


B Magnetic field [G] 




2x10^ 


n Particle density [cm^-^] 




5x 10" 



3.3. Particle distribution 

We assume an injection function Q{E) that is a power-law, 
of index 2.2, in the energy of the particles (Drury 1983), 

e(£) = eo£""e-^/^"'". (3) 

The normalization constant Qq for each kind of particle is ob- 
tained from, 

QoE-^EdE, (4) 






1-l-a 



and 



'^"1 + 1/fl 



= V 



QoE-"EdE. 



(5) 



Here a is the hadron-to-lepton energy density ratio. We con- 
sider a = 100 as observed in cosmic rays in the solar neigh- 
borhood (Ginzburg & Syrovatskii 1964), Ljei is the power in 
the form of relativistic particles and V is the acceleration vol- 
ume, i.e. the volume of the flux tube. The steady state particle 
distributions N(E) are obtained solving the transport equation 
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in the homogeneous approximation. 
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where fesc = fconv- This equation has an exact analytical solu- 
tion (see Ginzburg & Syrovatskii 1964): 



N(E): 



//' dE'Q(E') 
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with 



t(E,E')= / dE" 
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Figures |5] and |4] show the steady state particle distribution 
for protons and electrons. These distributions are valid only as 
long as we consider intervals much shorter than the timescale 
of the flare. 



Proton distribution 
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Log(E/eV) 

Fig. 4. — Steady state particle distribution for protons. 




Fig. 5. — Steady state particle distribution for electrons. 



3.4. Luminosity 

In order to estimate the luminosity due to electrons we com- 
puted synchrotron, relativistic Bremsstrahlung, synchrotron 
self-Compton (SSC) radiation, and IC up scattering of exter- 
nal seed photons (X-ray and IR radiation fields). We calculate 
for protons synchrotron and 7-ray emission from tt" decay in 
pp inelastic collisions. We consider also the synchrotron con- 
tribution from secondary pairs injected by charged pion decay 
(e.g. Orellana et al. 2007). The maximum particle density is 



^ 5x 10 cm . It is localized in the accretion columns. We 
consider a cylindrical portion of the accretion column of ra- 
dius tac ^ 10'° cm and height ~^ 0.1 r^c (e-g- Orlando et al. 
2010). This volume is where most of the relativistic interac- 
tions with matter take place^. For more details on high-energy 
processes see Vila & Aharonian (2009) and Vila (2010). 

3.4.1. Internal absorption 

Attenuation of the radiation produced by photon annihila- 
tion is expected in T Tauri stars. The opacity produced by a 
photon field with density nph(e) and photon energy e is 



t{E^) = 



1 
xdudedl. 



(1 - M)cr^^(£'^, e, M)Mph(e) 



(9) 



Here u = cos -d, ■& is the angle between the momenta of the 
colliding photons, / is the photon path, and a^^{E^,e,u) is 
the cross section for photon annihilation (Gould & Schreder 
1967). 

The absorbing photon fields are those generated within the 
system and the strong blackbody radiation field from the sur- 
roundings: the disk (IR), the star, and the X-rays from the ac- 
creting plasma. To estimate the disk infrared photon field we 
adopt an internal radius for the disk ofRu ^ 120 AU (see Fig. 
[B, and a temperature of T - 30 K (Dutrey et al. 1994). The 
X-ray temperature is taken to be T ~ 10^ K (e.g. Feigelson & 
Montmerle 1999). Regarding the star, we consider a tempera- 
ture of r^ ^ 4x 10^ K, and a radius R^, ^2 Rq. The opacity r 
(attenuation e""^) as a function of E^ is shown in Fig. |6l The 
absorption is almost complete above 100 GeV 
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Fig. 6. — Opacity curve as a function of energy E. 



4. APPLICATION TO POSSIBLE 7-RAY EMITTING T TAURI STARS 

After performing several catalog cross-correlations, a new 
possible association with nearby T Tauri stars inside the well- 
known p Ophiuchi star forming region clearly emerges for 
the Fermi source IFGL J 1625. 8 -2429c. Indeed, inside the 
95% confidence error ellipse of this Fermi source we find four 
T Tauri stars: 2MASS J 16260 160-2429449 (Casanova et al. 
1995), 2MASS J16253958-2426349 (Wilking et al. 1989), 

* The accretion disk is truncated at the edge of the magnetosphere. The 
exterior disk material has no major impact in the calculations since the inter- 
acting matter is in the accreting columns that penetrate the magnetosphere. 
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Right ascension 

Fig. 7. — GLIMPSE infrared image in the 8.0 ^m band sliowing the con- 
tents of tlie IFGL J1625.8— 2429c error circle towards tlie p Opliiuchi cloud. 
Several T Tauri stars are consistent with the Fermi 7-ray source position. 
They are labelled from 1 to 4 in decreasing order or right ascension. Axis 
coordinates are of equatorial J2000.0-type. 
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Fig. 9. — Radio mosaic of the IFGL J1625.8-2429c en'or circle as ob- 
served with the VLA at the 6 cm wavelength. Numbers indicate the T Tauri 
stars consistent with this Fermi source in decreasing order of right ascension. 
Only the T Tauri stai' # 3 is detected in the radio. 
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Fig. 8. — Composite X-ray image of the IFGL J1625.8— 2429c error circle 
obtained with the Chandra ACIS camera in the energy range 0.1-10 keV 
(Dataset identifier: ADS/Sa.CXO#obs/00618). Numbers indicate the T Tauri 
stars consistent with this Fermi source in decreasing order of right ascension. 
All of these stars are X-ray emitters. 

JCMTSF J162556.8-243014 (Di Francesco et al. 2008), and 
2MASS J16255752-2430317(Grasdalenetal. 1973). Then- 
observational properties are summarized in Table |3] In Figs. 
I2l|8]and|9]we show this field (/ =353?0, b =17?0) 

as observed in the infrared. X-ray and radio wavelengths. 
These images have been produced using public data retrieved 
(and calibrated when necessary) from the 5/?/teer-GLIMPSE, 
Chandra and NRAO Very Large Array (VLA) archives, re- 
spectively. As discussed below, we tentatively suggest that 
this Fermi source might be the result of the emission of at 
least these four Tauri stars that lay inside the location error 
box of IFGL J1625.8 -2429c. 

In order to estimate the a priori probability of a pure chance 
association we have implemented Monte Carlo simulations 
of computer-generated Fermi sources following the approach 
developed by Romero et al. (1999) for unidentified EGRET 
sources. After 10"* simulations of artificial Fermi populations 
we find 47 coincidences at 1 -degree binning and 4 at 2-degree 



1 



TABLE 2 

Fermi INTEGRATED FLUXES AND CALCULATED FLUXES. 

Energy range Flux ± error [cm"'^ s~' ] Computed flux 



100MeV-300MeV 
300MeV-l GeV 
1 GeV-3 GeV 
3 GeV- 10 GeV 
10 GeV- 100 GeV 



1.9x10"^ ± l.Ox lO"** 
2.2x10-'* ±0.6x10-** 
4x10-' ±0.7x10-' 
8.5x10-'° ±2.2x10 
2.1x10-'° ± l.Ox 10- 



10 



6.30x10-** 

1.89x10-** 

4.82x10 

1.21x10" 

3.76x10- 



1-9 



binning , indicating a probability of chance association of ^ 
lO""*. These results do not change with larger samples (we run 
up to 10^ simulations). 

Let F(Ei,E2) ± A[F{Ei,E2)] be the integrated Fermi flux 
in the energy range [Ei,E2] and its error To reproduce the 
observed fluxes, we consider in first approximation that the 
four T Tauri stars emit the same 7-ray luminosity. Then, we 
can compute the integrated flux from: 



F = 4 



WE) 
47rt/2£2 



AE. 



(10) 



Here Lj{E) is the total 7-ray luminosity produced by an indi- 
vidual T Tauri star according to our model (E > 20 MeV) and 
d is the distance to p Ophiuchi cloud, ^ 120 pc (Loinard et al. 
2008). The integrated Fermi fluxes (Abdo et al. 2010) and the 
calculated fluxes in five energy bands are shown in Table |2] 
Figure [To] shows the luminosity obtained with the model and 
the upper bound given by Chandra and Fermi data. Radio data 
from VLA are also shown, but this radiation is surely a com- 
bination of non-thermal and thermal emission and hence must 
be considered just as an upper value to constrain the model. 

5. DISCUSSION 

We see that for a reasonable set of physical parameters of 
T Tauri stars, a relatively weak 7-ray source can be produced. 
The intrinsic luminosity is not high {L^ '--^ 10^' erg s"' at ^ 1 
GeV), but since the stars are very nearby the flux is signif- 
icant. If such an association is confirmed, the T Tauri stars 
in p Ophiuchi would be the nearest 7-ray sources to the solar 
system detected so far. A caveat is in order, however, con- 
sidering the potential effects of the galactic cosmic rays with 
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TABLE 3 
Properties of T Tauri stars proposed as 1 FGL J 1 625 . 



-2429c counterpart. 



Star* 



Name 



Mag. 


X-ray counterpart 

cxo 


0.5-7 keV flux 
(erg s"' cm"^) 


6-cm flux density 
(mJy) 


R=\1.2 
V = 16.4 
i?=16.8 

y = i9.i 


J162601.6-242945 
J162557.5-243031 
J162556.0-243014 

J162539.5-242634 


2.30 X 10"'^ 
4.52 X 10-13 
4.95 X 10-'3 

4.90 X 10-1'* 


0.30 ±0.03 



1 2MASS J16260160-2429449 

2 2MASSJ16255752-2430317 

3 JCMTSFJ162556.8-243014/ 
2MASS J16255609-2430148 

4 2MASS J16253958-2426349 
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Fig. 10. — Computed non-thermal luminosity and Ferai! upper bounds for the four T Tauri stars, assuming a distance of 120 pc. The spectral energy distribution 
is connected by photon absorption. Model parameters as in Table 1 . 



the molecular material in p Ophiuchi. This emission has been 
subtracted from the Fermi 7-ray data (Abdo et al. 2010), but 
there could be more matter than detected so far. 

In the case that matter is underestimated in p Ophiuchi, 
we calculate the required density excess to explain the ori- 
gin of the high-energy emission from the source IFGL 
J1625.8-2429c. With this aim we assume that the whole 7- 
ray source is the result of passive "illumination" of cloud ma- 
terial by cosmic rays (e.g. Aharonian & Atoyan 1996). We 
adopt a standard cosmic-ray spectrum in p Ophiuchi. Proton- 
proton inelastic collisions with the ambient medium are the 
main radiative process contributing to the high-energy elec- 
tromagnetic spectrum. The cosmic -ray proton spectrum has 
a power law of index a = 2.7 and the enhancement parame- 
ter is taken to be 1 (i.e. absence of local acceleration). The 



normalization constant is obtained from the energy density 
of cosmic rays in the solar neighborhood, i.e. ujcr '^ 1 eV 
cm"^ (e.g. Ginzburg & Syrovatskii 1964). With these figures 
and the Fermi flux we obtain an average density of («) '^10'* 
cm"^, one order of magnitude greater than current estimates^ 
(e.g. Crutcher 1991). 

If the effects of the galactic cosmic rays with the interstel- 
lar material are responsible for the Fermi source, then IFGL 
J1625.8-2429c would be the nearest passive 7-ray source de- 
tected outside the solar system. 

6. CONCLUSIONS 



^ We remark that the cosmic-ray energy density in p Ophiuchi is not known 
with any accuracy, a fact that has impact in the estimates. However, a differ- 
ence of one order of magnitude would be, at least, peculiar. 



Are T Tauri stars gamma-ray emitters? 



We have found that under some assumptions T Tauri stars 
might be responsible for some nearby Fermi sources. We have 
presented a simplified model for the high-energy emission of 
this type of stars, that agrees with the available multiwave- 
length observations. T Tauri stars might be a new class of 
galactic 7-ray sources in the Galatic plane. Based on this new 
scenario, IFGL J1625.8-2429c is the first candidate for col- 
lective 7-ray emission from low-mass protostars. If the asso- 
ciation with the p Ophiuchi cloud is confirmed, it would be 
the closest 7-ray source to the Solar System. This statement 
still holds even in the alternative case where the detected 7- 
rays were simply due to cosmic rays interacting with the cloud 
ambient gas. 

More complex models can be developed in the future, us- 
ing new and deeper observations. The energy where the cutoff 
lays might be established in the future by new Cherenkov in- 
struments like CTA. 
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